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Abstract 

By reanalysing transverse momentum dependence in the perturbative calcula- 
tion of pion form factor an improved expression of pion form factor which takes 
into account the transverse momentum dependence in hard scattering ampli- 
tude and intrinsic transverse momentum dependence associated with pion wave 
functions is given to leading order, which is available for momentum transfers 
of the order of a few GeV as well as for Q — > oo. Our scheme can be extended 
to evaluate the contributions to the pion form factor beyond leading order. 
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1 Introduction 



Exclusive processes in perturbative quantum chromodynamics (PQCD) were first 
studied by Brodsky and Lepage J[| many years ago. Recently, the studies of the ex- 
clusive processes at experimentally accessible momentum transfers in the framework 
of PQCD have received much interest. The statement |2|, [| that the applicability 
of PQCD to the exclusive processes at experimentally accessible momentum trans- 
fers is questionable has been challenged [4-11]. Huang and Shen || pointed out that 
the applicability of PQCD to form factor is questionable only as momentum transfer 
Q 2 < 4 GeV 2 through the study on the pion form factor by reanalyzing the contribu- 
tions from end-point regions. Botts and Sterman || proposed a formulation in which 
both Sudakov and nonleading logarithmic corrections to independent (Landshoff) 
scatterings |12| of valence quarks are organized systematically. Li and Sterman || 
gave out a modified expression for the pion form factor by taking into account the 
customarily neglected partonic transverse momentum as well as Sudakov correction. 
They reached a similar conclusion as [[|: PQCD begins to be self-consistent at about 
Q ~ 20Aqcd- Jakob and Kroll || demonstrated that for momentum transfers of the 
order of a few GeV the intrinsic transverse momentum dependence of wave function 
leads to a substantial suppression of the perturbative contributions, which should be 
considered besides Sudakov suppression. The base of the most previous discussions 
[5-11] is the formalism in Ref. || which is suitable for studying the large Q region 
since it sets b — > in the integral of the wave function in respecting the intrinsic 
transverse momentum dependence of wave function. In this paper we re-analyse the 
PQCD calculation for the pion form factor at momentum transfers of the order of 
a few GeV. We will give out an improved expression for the pion form factor which 
takes into account the transverse momentum dependence in the gluon propagator 
as well as in the fermion propagator in the hard scattering amplitude T# and the 
intrinsic transverse momentum dependence associated with the pion wave function. 
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This expression is available for momentum transfers of the order of a few GeV as 
well as for Q — > oo since it dose not make the approximation b — > 0. The formalism 
of Ref. || is just an approximate expression in respecting the intrinsic transverse 
momentum dependence. This approximation brings sizeable effect on the numeri- 
cal prediction for the pion form factor in the momentum transfer Q ~ a few GeV 
region. It is also found that the transverse momentum dependence of the fermion 
propagator in T H leads to a mild reduction of the prediction for the pion form factor 
in the same momentum transfer region. The remainder of the paper is organized as 
follows. Sect. 2 reviews Ll-Sterman's formalism for the pion form factor. In Sect. 
3, we discuss the transverse momentum dependence in pion wave function beside 
the one in the hard scattering amplitude T# in perturbative calculation. In Sect. 4, 
we do our numerical calculations. Finally, in Sect. 5, we give a summary. 



2 Brief review of Li-Sterman's formalism 

Taking into account the transverse momenta k^ that flow from the wave functions 
through the hard scattering leads to a factorization form with two wave functions 
if>(x{, kyj corresponding to the external pions, combined with a new hard-scattering 
function Tjj(xx, x 2 , Q, ky 1( k^J, which depends in general on transverse as well as 
longitudinal momenta ||, 

FAQ 2 ) = I dxidx 2 f „ kT ! d ^ if) (xi , k Tl ) 

x T H (xi,x 2 ,Q, k Tl ,k T2 )V(^2,kT 2 ). (1) 

In this form, both soft and collinear logarithmic enhancements are factorized into 
the functions if). At the lowest order, T# reads 

T H {x u x 2 ,QM,M.) - {xiQ 2 + k 2 TiHxiX2Q 2 + ikTi+kT2n (2) 



2 A frame with qr = has been adopted to obtain this expression of Th in @ . The expressions 
of Tu with different momentum assignment will be analysed in detail |fl3[ . 



Neglecting the transverse momentum dependence in the fermion propagator, Eq. 
(H) becomes 

167rCVa s (/x) 



T H (x 1 ,x 2 ,Q,'k Tl ,kT 2 ) 



(3) 



X1X2Q 2 + (k Tl + k T J 2 ' 
The next step is to re-express Eq. (P in terms of the Fourier transformation 

variables in the transverse configuration space ||. Observing that Th in Eq. (§) 

depends on only a single combination of the transverse momenta (k^ + k^ 2 ), the 

Fourier transformation of Eq. (TO) involves only a single integral of Fourier transform 

variable 6, 



F^Q 2 







1 dh 

dx 1 dx 27 -—-^(p(xi, b, n)T H (x 1 ,x 2 , Q, b, fi)(p{x 2 , b, /i). (4) 

The wave functions in 6-space, (p(x, b, /i) take into account an infinite summation of 

higher-order effects associated with the elastic scattering of valence partons, which 

gives out Sudakov suppression to the large-6 and small-x regions. 

The asymptotic behavior of <p(x, b, /i) at large Q 2 has been obtained in Ref. |5j 

w dfi 



ip(x, b, fi) = exp 



x 



1 

X,- 



(5) 



[i 1 dfx 

-s{x, b, Q) - s(l - x, b, Q) - 2 / —j g (g(fl)) 

Ji/b fi 

where 7 9 = —a s /7c is the quark anomalous dimension in the axial gauge. s(^,b,Q) 
is Sudakov exponential factor, which reads |5], |TT| 



s(£,b,Q) 



AW) 
28, 
A^3 2 



q In 
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m 



ln(-26) + 1 
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ln(-2g) 



2/?i 
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8^ 3 



ln 2 (2g) - In 2 (-26) 



(6) 



where 



ln[£Q/(V2A)], S = ln(6A), 
33 - 2n f 153 - 19n f 
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a(2) 67 1 2 10 8 ^ , .1 „, 

= 7T 2 rif + -Bi h(-e 7 

9 3 27 t 3^ V 2 ; 



(7) 



rif is the number of quark flavors. 7 is the Euler constant. The factor e~ s ^' b '^ 
which induces an enhancement in small-6 regions has been set to unity whenever 
£ < y2/(bQ). The coefficients in Eq. (El) are different from Ref. since there are 
some algebraically mistakes in the previous expression of s(£, b, Q), which has been 
pointed in Ref. |IT| . (x, |) in Eq. ([|) is a "soft" wave function calculated with 



gluons of transverse momentum fcr < 1/6, 



Neglecting the 6-dependence of the function 0(x, 1/6) and performing Fourier 
transformation for Eq. (^), the expression Eq. (|TJ) becomes || 

/*1 /*oo 

F W (Q 2 ) = IQirCp / dxidx 2 / 6 dba s {t)KQ{^/x~[X2Qb) 
Jo Jo 

x exp(-S , (xi,x 2 ,Q,6,t))0(a;i)0(a;2), (9) 



where is defined as 



d 2 k 



m= i ^^(x,k T ), (10) 



_K"o is the modified Bessel function of order zero, and 
S(x l ,x 2 ,Q,b,t) 



2 t 

^ (s(xi, b, Q) + s(l - x i; b, Q)) - — In — * 
»=i Pi ~b 



3 /cy-dependence in perturbative calculation 

Following Ref. 0, we find that 4>(x, 1/6) in Eq. (H) should be replaced by the 
amplitude 

(p m^ b )=f ^ e ik ^(x,k T ). (12) 

Jk T <l/b l07T d 

Through the requirement that e ikr ' b ~ 1, as Q — > 00 (6 — > 0). Ref. || proposed 
ip^ b >(x, b) can be expressed approximately by (p(x,l/b). For simplicity, Ref. || 
neglected the 6-dependence of the function Lp^^ b \x, b), namely neglected the intrinsic 



transverse momentum dependence of wave function. Then ipw b '(x, b) can be pulled 
out of the fe-integral. Jakob and Kroll || improved this approximation and proposed 
to take into account the intrinsic transverse momentum dependence by a function 
i/j(x, b) which is the Fourier transformation of wave function, 

iP(x,b)=[ ^ e ^V(*,k T ). (13) 
Jk T <oo Ibn 6 

It can be seen that ip(x, b) is just an approximate expression for ,b), i.e. 

ip(x,b) = (p(°°\x,b). When b ^ (namely Q — > oo) they are consistent with each 
other. When Q is the order of a few GeV, the difference may be sizable. The 
difference between cp^' b '(x, b) and ip(x, b) can be expressed by a function D(x, b), 

poo rf^lfrp 

D(x,b) = j i/bT ^j^(x,k T ), (14) 

which increases as b becomes large. Eq. ( |TB| ) enlarges the upper limit of the integral 
Eq. ( |T2D from 1/b to oo, which corresponds to evaluate the contributions from 
the perturbative tail of wave function once again. Sudakov form factor provides 
much more suppression for large Q than for small Q. Thus substituting ip(x,b) for 
(pV-' b >(x,b) does not effect the pion form factor for large Q. As Q is the order of a 
few GeV one should investigate the effects of this substitution. 

The Fourier-transformed hard-scattering amplitude from Eq. (|j) reads 

T H (xi,x 2 , Q, b, fi) = l&Ka 8 (fi)C F K (^/xix 2 Qb) . (15) 



The renormalization group applied to T# gives 



-4 f d 4l q {gm 



xT H ( Xl ,x 2 ,Q,b,t). (16) 



T H (x 1 ,x 2 ,Q,b,fi) = exp 



The variable t was taken as t = max ly/xiXzQ, 1/b) and a cut-off on the running 
coupling constant (a s < 0.7) for large-6 region was made to guarantee PQCD to 
be self-consistent in Ref. 0. However, in the regions of small XiX 2 Q 2 and large b, 
the nonperturbative contributions maybe important. For example, the multi-gluon 



exchange can occur between quark and antiquark and the transverse momentum 
intrinsic to the bound state wave-functions flows through all the propagators. Ref. 
14| suggested a frozen ot s to take into account these effects. Instead of the cut-off 



of a s in the 6-space, Ref. |TT[ suggested that the coupling constant is frozen at 



b ~ l/(y (k^)) and the variable t is taken as 

t = max (y/x\X 2 Q, 1/&f) , (17) 

where 



if 1/6 > J (k|) 

~ V ~ (18) 



1/^04) if i/6< 

and -y/ (k^) is the average transverse momentum of the pion. In this way, the per- 
turbative contributions to the pion form factor can be calculated from the present 
energy with a reasonable a s . It should be emphasized that the average transverse 



momentum J (k^) is determined definitely by the hadronic wave function. 
Combining Eqs. (|), (g), (|TJ), (|TJ) and (|TJ), we have 

Ftt(Q 2 ) = IQttCf / dx\dx 2 \ b dba s (t)K (^xix 2 Qb) 
Jo Jo 

x exp (S(x l ,x 2 , Q, b, t)) V {1/b \x 1 , b) V ^ b \x 2 , b), (19) 



where 

S(x 1} x 2j Q,b, t) 



2 t 

(s(xi, b, Q) + s(l - Xi, b, Q)) - — In 
i=i Pi -6 



(20) 



Eq. (|T^) is an improved expression for the pion form factor. From our formalism, it 
can be found that: (a) neglecting the transverse momentum dependence associated 
with the wave function, Eq. (^) becomes the expression in Ref. |5| (See Eq. (^)); 
(b) approximating (p^-^ b '(x,b) with ?/>(x,&), Eq. flT9"|) becomes the expression in Ref. 

I, 

F-k(Q 2 ) = IQttCf / dx\dx 2 I b dba s (t)K (y / xiX2Qb) 
Jo Jo 

x exp (-S(xi,x 2 ,Q,b, t))ip(x 1} b)ijj(x 2 ,b), (21) 



where ip(x, b) is the Fourier transformation of the wave function (Eq. fll3D ). 

The /c^-dependence in the fermion propagator contributes to T H a factor Zl Q2^ k ^ 
which involves only a single transverse momentum corresponding to one in the ex- 
ternal pion (see Eqs. (EJ) and (|3])) and this factor makes the Fourier transformation 
for Th involve multi-6-integrals0. Although the /^-dependence in this factor is 
linear rather than quadratic in the x's, it may bring some effects at the end-point 
region x\ ~ 0. We can consider this factor and keep the simplification of the Fourier 
transformation for Tjj at the same time by combining it with the corresponding 
wave function to define a new "wave function" , 

xiQ 2 

ip(x u k Tl ) = ————$(x u k Tl ). (22) 

Substituting tp for ip in Eq. (|1|), the transverse momentum dependence of the fermion 
propagator in T H can be taken into account easily in the perturbative calculations 
for the pion form factor. 



4 Numerical calculations 



We adopt two models of the pion wave function: (a) the BHL wave function |H | 



ip {a) (x,k T ) = A exp 



k£ + m 2 



l(3 2 x(l - x) 



(23) 



where |TJ A = 32 GeV" 1 , (3 = 0.385 GeV, m = 289 MeV and yj (k 2 T ) = 356 MeV; 
(b) the CZ-like wave function [16-18] 

kl + m 2 



4> {b \x,k T ) = A (1 - 2xf exp 



^ 2 x(l - x) 



(24) 



where pf A = 136 GeV" 1 , p = 0.455 GeV, m = 342 MeV and ^(k 2 .) = 343 MeV. 

In order to discuss the effect of transverse momentum dependence associated 
with wave function, first, we neglect the /^-dependence of the fermion propagator 
in Th. Fig. 1 compares the behaviors of functions ip^-' b ^(x, b), ip( x , b) and (f)(x, 1/6). 
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All of them suppress the the contributions from large-6 region, but the suppression 
behaviors are different in quantity. In the b ~ region, <f)(x, 1/6) and tp(x,b) are 
good approximation for ip^' b '(x, 6), while in the large-6 region that is questionable. 
Sudakov form factor suppresses the contributions from large-6 region in a mild way 
for small momentum transfer Q than for large Q. Thus it is not a good approxi- 
mation with tp(x, 1/b) to respect the intrinsic A/p-dependence of the wave functions 
for small Q, since much more contributions come from large-6 region. Numerical 
evaluations of F n through Eqs. (||), (|19"D and (pT|), using the BHL and CZ-like 
wave functions, which are plotted in Figs. 2 and 3 confirm the observation made in 
Fig. 1. Comparing with the original PQCD prediction fl[ of the pion form factor 
(namely neglecting the /c^-dependence in T H as well as in wave functions; the dotted 
line), Sudakov form factor gives a suppression effect to F^, which was first pointed 
out by Brodsky and Lepage |TJ. Comparing with the result obtained from Eq.(|^) 
(namely neglecting the intrinsic fey-dependence of wave function; the dash-dotted 
line), b) (the solid line) suppresses the contribution from PQCD by about 

50%, while ip(x, b) (the dashed line) suppresses by about 30% at Q = 2GeV for the 
BHL wave function. The corresponding quantities are 55% and 35% respectively in 
the case of the CZ-like wave function. The intrinsic k^- dependence of the wave func- 
tion provides additional substantial suppression for F n besides Sudakov form factor. 
ip l l b (x,b) suppresses the contribution from PQCD more strongly than ?p(x,b) does. 
The numerical predictions of F n obtained from b) and ijj(x,b) are different in 

the momentum transfer Q ~ a few GeV region. 

We evaluate the effect of the /cr-dependence in the fermion propagator also. It 
leads to a small reduction of the prediction for F n by about 10% for Q ~ 3 GeV, 
which is coincides with Ref. 0. 
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5 Summary 



In this paper we re-analyze the transverse momentum dependence in the perturba- 
tive calculation of pion form factor at momentum transfers of the order of a few 
GeV. We give out an improved expression for the pion form factor which takes 
into account gluon propagator as well as fermion propagator transverse momentum 
dependence in the hard scattering amplitude and intrinsic transverse momentum 
dependence associated with pion wave functions. It is found that the previous ap- 
proach is just an approximate expression in respecting the transverse momentum 
dependence associated with wave functions. This approximation brings sizeable ef- 
fect on the numerical predictions for the pion form factor in the momentum transfer 
Q ~ a few GeV region. It is also found that the transverse momentum dependence 
of fermion propagator in Th leads to a mild reduction of the prediction for the pion 
form factor in the same momentum transfer region. 

We would like point out one more times that our formalism is available for 
momentum transfer of the order of a few GeV as well as for Q — > oo, and our scheme 
can be extended to evaluate the higher order and higher helicity contributions for 
the pion form factor. The more studies on the /^-dependence in the hard scattering 
amplitude Th are in proceeding. 
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Figure Captions 

Fig. 1. The functions (p^ b \x,b) (solid line), ip(x, 1/6) (dashed line) and 4>(x,b) 
(dotted line) which are adopted to respect the intrinsic transverse momentum 
dependence. 

x = 0.5, Q 2 = 4 GeV 2 , for the BHL wave function. 

Fig. 2. The pion form factor calculated with the BHL wave function. The solid 
and dashed lines are for ip^^ b \x, b) and ip(x, 1/6) respectively. The dash-dotted 
line is obtained by neglecting the intrinsic transverse momentum dependence 
in wave functions. The dotted line is evaluated by neglecting fc T - dependence 
in both Th and wave functions. 

Fig. 3. Similar to Fig. 2. The pion form factor calculated with the CZ-like wave 
function. 
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